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Asymptotic expansion and central limit 
theorem for quadratic variations of 
Gaussian processes 
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Cohen, Guyon, Perrin and Pontier have given assumptions under which the second-order 
quadratic variations of a Gaussian process converge almost surely to a deterministic limit. In 
this paper we present two new convergence results about these variations: the first is a determin- 
istic asymptotic expansion; the second is a central limit theorem. Next we apply these results to 
identify two-parameter fractional Brownian motion and anisotropic fractional Brownian motion. 
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Introduction 



In this paper we consider second-order quadratic variations of a Gaussian process X. 
From Cohen et al. [10] we know that it converges to a deterministic limit under convenient 
conditions on the covariance function of the process. First we sharpen this result: we 
show that if an asymptotic expansion of the covariance function is known, we get an 
asymptotic expansion of the second-order quadratic variation. Next we establish a central 
limit theorem related to the previous result. We apply these results to two-parameter 
fractional Brownian motion, which is a generalization of fractional Brownian motion 
that has non-stationary increments, and to anisotropic fractional Brownian field, which 
is a multidimensional anisotropic generalization of fractional Brownian motion. 

In the first section, we state the notation. In the second section, we prove the main 
theorems about the second-order quadratic variation. In the third section, we study the 
case of two-parameter fractional Brownian motion. In the fourth section, we consider 
anisotropic fractional Brownian motion. 
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1. Notation 
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Let X = {X t ;t £ [0, 1]} be a Gaussian process. We denote by 1 i— > M t its mean function 
and by (s,i) i— * R(s,t) its covariance function. 
We define the second-order increments of R as 

S^R(s, t) = R(s + h,t) + R(s -h,t)- 2R(s, t), 
5%R(s, t) = R(s, t + h)+ R(s, t-h)- 2R(s, t). 

For fractional processes (i.e., processes whose properties are close to those of fractional 
Brownian motion), we use the second-order quadratic variation 

n-l 

Vn(X) = y^[X( fc+ i)/ n +X( k _ 1 y n - 2X k / n ] , (1) 
fc=l 

because the standard quadratic variation does not satisfy a central limit theorem in 
general. 

To be sure that V n (X) converges almost surely to a deterministic limit, we need to 
normalize this quantity. A result of the form 

lim n 1 ~ 1 V n (X) = [ g{t)dt a.s. (2) 

n-»+oo J Q 

is expected, where 7 is related to the regularity of the paths of X, and g is related 
to the non-diffcrcntiability of R on the diagonal {s = t} and is called the singularity 
function of the process. In this paper, we consider a class of processes for which a more 
general normalization is needed. Moreover, we present a better result because we give an 
asymptotic expansion of the left-hand side of (2). 

We will say that a Borel function ip : ]0, a[— ► M (a > 0) is regularly varying with index 
(3 € M if ip(h) = h^L(h), where L is a slowly varying function 

VA > lim ^ = 1. 
.t^o+ L(x) 

Let deW. Standard fractional Brownian motion (FBM) B H = {B?;te R d }, with Hurst 
index H s]0, 1[, is the unique continuous centered Gaussian process, which has the co- 
variance function 

Vs,teR d Cov(Bf, J Bf) = |(|s| 2H + \t\ 2H -\s-t\ 2H ), (3) 

where | • | denotes the Euclidean norm. 

In next section, we use the notation (we drop the superscript index n wherever it is 
possible) 

AX fe ™ =X( k+1 y n + X( k _ 1 y n -2X k / n , k = l,...,n — 1, 
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and 

d&> = E(AX^AX^), j, k = 1, . . .,» - 1. (4) 

2. The results 

In this section we sharpen (2). First, we prove a deterministic asymptotic expansion of 
V n (X) under certain conditions on the covariance function. Second, we prove a central 
limit theorem. 

Examples of the application of Theorem 2.1 with a non-trivial slowly varying function 
L(h) can be found in Section 4.2. 



2.1. Asymptotic expansion 



Theorem 1. Assume that X satisfies the following statements: 

1. t\—*Mt = KX t has a bounded first derivative in [0,1]. 

2. The covariance function R has the following properties: 

(a) R is continuous in [0, l] 2 . 

(b) The derivative g f2 exists and is continuous in ]0, l] 2 \ {s = t} . There exist a 
constant C > 0, a realj s]0, 2[ and a positive slowly varying function L:]0, 1[— > 
]0, +oo[ such that 



Vs,te}0,l} 2 \{s = t.} 



d A R 



ds 2 dt 2 



(s,t) 



<C- 



.L(\s-t\ 



t |2+ 7 - 



(5) 



(c) There exist q+l functions (q € N) ga,gi, ■ ■ ■ ,g q from ]0, 1[ to R, q real numbers 
< v\ < ■ ■ ■ < i/ q and a function 4> : ]0, 1[— > ]0, +oo[ such that: 

(i) if q •_ 1, then V0 < i < g — 1, is Lipschitz on ]0, 1[; 

(ii) <7 g is bounded on ]0, 1[; 

(iii) we ftave 



sup 

/i<t<i-ft, 



h 2 -~<L(K) 



-g Q (t)-J2^(t)Hhr 



/l = 0+ o(</>(/^), (6) 



where if q = 0, i/ien Xh=i ffiW^W* = a7 ^ 4>{h) Uq = 1; efee if q ^ 0, 
then lim/ l _ > o+ </>(/*) = 0. 
3. Ifq^O, we assume that 

logn 



lim 



0. 



n— >+oo ncj)(l/n) v i 
4. If Xis not centered, we make the additional assumption 

1 



lim ; ; — ; ; — ; — ; = 0, 

n^+oo niL(l/n)(j){l/nyi 



(7) 



(8) 
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where if q = 0, then <f>(l/n) Vq = 1 . 
Then, for all t € [0, 1], we have almost surely 



WR, 



■Vn(X) 



g (x)dx + ^2( 9i(x)dx Ul - j +o 



(9) 



Remarks, (i) If the assumption (6) is fulfilled for g*, then it is fulfilled for all q£ 
{0,1, ... ,q*} too with the truncated sequences (<?i)o<i<g and (vi)o<i< q . The maximal 
value of q is given by the assumption (8), which yields an upper bound for the value of 
v q . 

(ii) Assumption 2 in Theorem 1 implies that the functions gi,0 < i < q, arc continuous 
and bounded on ]0, 1[, and so they are Riemann integrable on this interval. 

(iii) In the case 7 > 1, the assumption (8) is a consequence of the assumption (7) and 
of Karamata's representation of positive slowly varying functions (see Bingham, Goldic 
and Teugels [7], Theorem 1.3.1). 



Proof of Theorem 1. We set v = and fix the convention that <f>(K) v ° = 1. Moreover, 
in the entire proof K denotes a positive constant whose value does not matter. First we 
assume that X is centered. 

We prove the following asymptotic expansion for the expectation of V n (X): 



'( - 



'i 5 



(10) 



We have 



and 



d jk = (Sl /n o5l /n R)(i,*) 
\n n J 



n-l 



EV n {X) = J2dkk- 
Moreover, the assumption (6) yields 



fe=i 



sup 

k=l n-l 



n-y- 2 L(l/n) 



Therefore, 



71 — t' + CX) 



(11) 



(12) 



(13) 



lim sup — 



— EV n (X)- V / 9i {x)dx4> 
In) ~—iJc\ 



L{l/n) 
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< lim sup — 

n~>+oo (p{l/n) u 1 



Tl-1 

E 



Hk 



n — * \ 11 



+ limsup 

n^+oo <p{\/n) v i 



1 / j.1 1 n-l 

Ely ft(*)d«--E 

i=0 V° k=l 



n J \n 

k 
n 



We get 



L x = lim sup— — — 

n->+oo (f>{l/n) v i 



E 



dkk 



-X>(- 



71/ \ II 



< lim sup , — sup 

„_> +00 0(l/n) l, 'f fe=i,..., n _i 



'A A 



n~t- 2 L{l/n) 



i=0 



Thus (13) implies that Li = 0. 
For L2, notice that 



1 



L 2 = lim sup 

iw+oo <p(l/n) l '9 



< lim sup — y 



-V 



i=0 



r (Ti-l)/?i 



n-l 



9 



:(x)dx E^l ~ 



k=l 



- lim sup 



n— *+oc 



<j>{l/n) v " 7 ( „_i)/„ 



54 (.T)d a; = L 2 1) +4 2) . 



The term is obviously equal to due to (7) and the fact that the functions gi are 
bounded. 

Moreover the assumption 2(c) (i) in Theorem 1 implies that there exists K > such 
that for all 0<i<q- 1, 



(n-l)/n 



n-l 



3 







n-l fkf 

^E 



1 J(k-X)/n 



(x)dr--^ 

A-=] 

0*0*0 - .9 



n-l /.fe/n 

dx<^V / 

k=1 J(k-l)/n 



— x ) dx < 



K 



Consequently, 



L^limsup 



71 — >-+OQ 



</>{i/ n y 
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< lira sup , ,„ , . — 



71 — > + CX> 



/n) 



n — > + do 



- lim sup 

9-1 



i=0 

p(n-l)/n ^ "-1 

/ ff g (a;)cia; V.g f; 

•'O n fe=l 



n-1 



(n-l)/n ^ ^ ^ 

S>i(x)d:r VV 



fe=l 



< 



1 



X lim sup ; — — — 

^ n^+o5 n^(l/n)"« 

+ lira sup 



n — - t do 



In) 
(n-l)/n 



g q (x) da; - - V 5<3 ( - 
n \n 

fc=l v 



where the first term of the right-hand side is equal to because of (7) and the assump- 
tion 2(c)(i) in Theorem 1, and the second term is equal to according to classical results 
on Ricmann sums. Therefore, = 0. This proves the asymptotic expansion (10). 
Next we prove that almost surely 

„l-7 



: {V n (X)-EV n (X)) 



n— > + oo 



■I 5 



L(l/n 

Application of Cochran's theorem to the Gaussian vector yields 



(14) 



1,1-7 



0(l/n)"«L(l/n) 



So there are n — 1 nonnegative real numbers (Mi n > • • • >Mn-i,n) an( i one i n ~ 1)~ 
dimensional Gaussian vector Y n , such that its components are independent Gaussian 
variables Af(0, 1) and 



,1-7 



n-1 



L(l/n) 



(15) 



3=1 



As in Bcgyn [5], Hanson and Wright's inequality (sec Hanson and Wright [12]) yields 
that, for all < e < 1, 



,1-7 



I ^^—\V n (X)-EV n {X)\>e)<2e SP (-Ke 



no\ — 

. n 



(16) 



So if we set 
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it follows from (7) that 

+00 



JHm,e„ = and gp^ ___ |F»(X) - EV n (X)| > e n j < +oo, 

and the Borel-Cantelli lemma yields (14). 

Now let us examine the case of non-centered X. Set M = {M t ;i S [0, 1]}. From as- 
sumption 1, 

1 ' -V n {M) n ^°° o( 



L(l/n) v ' \n'fL(l/n) 
and by adding (8) we obtain 

n— >+oo L\\/n) 

So, if we apply the theorem to the centered process X t = — E(X t ), using the arguments 
of Baxter [3], we obtain the result for X . □ 

In the sequel, we apply these results to the identification of some fractional models. We 
will obtain strongly consistent estimators that will be more interesting in practice if they 
are asymptotically normal. Therefore, we establish a central limit theorem for V n (X). 



2.2. Central limit theorem 

The integral 



j+l i-u r k+l 

du j 



P7 (j,k)= 1 du 1 dv j k dx J x i {v _ y) 2 +1 d y 



1 



with j — k > 2 is absolutely convergent when 7 < 2. Because it depends only on the 
difference j — k, we denote it p 7 (j — k). 

By considering I = j — k > 2 and < 7 < 2, we obtain the following equalities: If 7 7^ 1, 

(\l 2| 2 ' 7 - 4|Z - 1| 2 ^ + 6|Z| 2 ^ 4\l + 1| 2 ^ + \l + 2| 2 ' 7 ) 
* (0 = (7 -2)( 7 - 1)7(7 + D ' ^ 

if 7 = 1, 

Pl (l) = |(|Z - 2| log \l-2\- 4|Z - l|log|Z - 1| + 6|i| log|2| 

- 4|Z + 1| log 1 1 + 1| + \l + 2| log \l + 2|). (19) 
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Moreover, we notice that (17) yields the existence of a constant K > such that, for all 
I > 2, we have |p 7 (/)| <Kl~ 2 -t. For 7G]0,2[, set 



I Ay I 



(20) 



We may now prove a central limit theorem with additional assumptions. The preceding 
formulas will be useful to compute the asymptotic behavior of djk- 

Theorem 2. Assume that X is centered and satisfies the following statements: 

1. R is continuous in [0,1] . 

2. Let T = {0 <t < s i< 1} . We assume that the derivative d ® 2 g t 2 exists in ]0, l] 2 \ {s — 
t} , and that there exist a continuous function C:T^R ; a real 7 6]0,2[ and a 
positive slowly varying function L : ]0, 1] — > K such that 



V(a,*)er 



(s-t) 2 +~> d 4 R 
L(s-t) ds 2 dt 2 



(s,t) = C(s,t), 



(21) 



where T denotes the interior of T (i.e., T= {0 < t < s < 1} ). 
3. We assume that there exist q + I functions (q S N) go,9i, ■ ■ ■ ,9q from ]0, 1[ to R, q 
real numbers < v\ < ■ ■ ■ < v q and a function 4>:]0, 1[— *]0, +00 [ such that: 

(a) if q > 1, then V0 < i < q — 1, gi is Lipschitz on ]0, 1[; 

(b) g q is (1/2 + a q )-H6lderian on ]0, 1[ with < a q < 1/2; 

(c) there exists t G]0, 1[ such that go(t) ^ 0; 

(d) we have 



lim 



'( 



^ . sup 

/w0+ V/i \fe<t<l-fc 



(<Jf o^J?)(t,t) 



h 2 -tL(h) 



»>(*) -!>(*) w 



0. 



(22) 



where if q = 0, f/ien X)?=i 9i{t)4>{l~ l Y i — 0; where if q 7^ 0, t/ien lim/j^o+ (f>(h) = 0; 
(e) £/iere exists a bounded function g :]0,1[— >R such that 



lim sup 

h^0+ h<t<l-2h 



(S'{o5%R){t + h,t) 



h 2 ^L(h) 



= 0. 



Then we have 



L(l/n) 



go (a) da: / g i (x)da; 



1=1 



■AA(0,a 2 



where 



g (x) 2 dx + 4 / 5(x) 2 dx + 4||p 7 || 2 / C(x, x) 2 dx 



(23) 



(24) 



(25) 



720 



A. Begyn 



Remark, (i) In Theorem 1, there was no minimum value for the integer q, but this is 
not the case in the assumption (22): we must choose q large enough such that the bias 
is negligible with respect to the stochastic error. 

(ii) Assumption (22) yields that the functions g.j, < i < q, are continuous on ]0, 1[. 

Proof of Theorem 2. In all the proof, K denotes a positive constant whose value does 
not matter. To simplify notation, choose the convention Vq = and Vh <E ]0, 1[, (j)(h) va = 1. 
Set 

bn = J2 9i(?)te<l> - , T n = ^i——V n {X), T n =T n — E(T n ). (26) 
f^Jo \ n J L ( l / n ) 

We split the proof into three steps: in the first and second steps, we prove the convergence 
when n — > +oo of T n toward a centered Gaussian law with variance a 2 ; in the third step, 
we prove the conclusion of Theorem 2. 

Step 1. We note that V n (X) is the square of the Euclidean norm of a (n— l)-dimensional 
Gaussian vector whose components are 

Hence, by the classical Cochran theorem, we can find a n = n — 1 positive real numbers 
(Ai jn , . . . , Xa n ,n) and one a„-dimensional Gaussian vector Y ni such that its components 
are independent Gaussian variables A^(0, 1) and 

V n (X) = -£\ j!n (Y(») 2 , (27) 
j=i 

with the convention that the empty sum is equal to zero. 
We set 

S n (X) = V n (X)-EV n (X). 
We want to apply the Lindcberg central limit theorem to S n (X). We must verify that 

A;= max \ jin = o(y/VaxS n (X)). (28) 

l<j<a n 

We have 

1--Y 

Th x — 

A* < K — ; — -— max > |d,-t|. 
" ~ L(l/n) i<k<n-ij-^ 1 Jkl 

With the same methods as in Begyn [5] , we can check that 
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We have 



E[V n (X) 2 ] = E E[(AX,) 2 (AX fc ) 2 ]. 

j,k=l 

Because the vector (AXk)i< k <N-i is Gaussian, Isserlis formulas yield (see Isserlis [14]) 

E[(Al t ) 4 ]=3(E[(AX,) 2 ]) 2 

and, if j 7^ k, 



Therefore, 



E[(AX 3 ) 2 (AX k ) 2 } = E[(AX 3 ) 2 n(AX k ) 2 } + 2(E[AX 3 AX k 



VaxV n (X) = 2j24k+* E d % 

fe=l l<k<j<n-l 



(29) 



and, consequently, 

VaxS n (X) > 2 E d ; 
Moreover, the assumption (22) yields 



2 

fefc- 



fc=i 



hm sup 

ra-H-00 i<fe< n _i 



and, because 50 is bounded on ]0, 1[, 



(S[ /n o5l /n R){k/n,k/n) (k 

— go 



?iT- 2 L(l/n) 



= 



lim sup 

"- , + 00 KKn-l 



(5\ /n o8y n R){k/n,k/n) 
n~<- 2 L(l/n) 



9o 



With the same ideas as in the proof of (10), we can show that the previous limit yields 

L{l/n) z f— J ,/n 



Therefore, 



,3-27 



lim inf . . N9 Var S„ (X) > 2 / 5 2 (x) dx>0. 

ni/ny J 

Thus there exists K > such that 

0< - = <* 
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which yields (28). 

Consequently, the Lindeberg central limit theorem yields that when n — > +00, 



S n (X) 



x /Var5„(X) 

Equivalently, we have shown that when n—* +00, 

T n (£) 



AA(0,1). 



AA(0,1). 



VarT„ 



Step 2. Let us prove that 



where a 2 was defined in (25). 
We have seen (29): 



lira Var T n = a, 

n — >-f 00 



l<k<j<n-l 



k=l 



Step 2.1. Let us prove that 

„3-2 7 

lim 



n->+oo L(l/n) 2 



»1 +oc 
<j<n-l 1=3 



<k<j<n-l 
j-k>3 



If j 7^ 1, fc 7^ 1 and j — k > 3, then 

\n n J 



(j + l)/n pu 

du 

j/n 



dv 



u—(l/n) Jk/ 



(k+l)/n fx 

dx 



d 4 n 



x-(l/n) ds 2 dt 2 



(v,y) dy 



(j + l)/n f-u 



du 



dv 



(fc+l)/n rx 

dx 



u-(l/n) Jk/n J x-(l/n) i v - V) 2+1 



C(v,y) 



(30) 



(31) 



(32) 



■L(v-y)dy. 



We set 
e„ = sup 



C(v,y)-C[ J -,- 



„ 7-1 7 + 1 ft-1 fc + 1 

3 < j — k<n-2,- < v < , < y < 

n n n n 



Because C is uniformly continuous on the compact set T, 



lim e„ = 0. 

n — >+oo 
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Moreover, we set 
r n {j,k) 



(j + l)/n f-u 

du 

j/n 



dv 



(fc+l)/n fx 

dx 



,-r-a 



j+i r u 



k+l 

dv I dx 



L(v-y) 

x _ ( i /n) («-J/) 2 +t 
* L((v-y)/n) 



dy 



x-i («-y) 2+ ^ 



and because it depends only on j — k, wc denote it r n (j — k). 
We have 



ijk 



C 



L * 

n' n 



(j+l)/n r « 
dlt 

j/n 

< e n r n (j - k). 



dv 



(k-\-l)/n c-x 

dx 



C(v,y)-C(j/n,k/n) 



u-(l/n) Jk/n Jx-(l/n) ( v ~ 2/) 2+7 



L{v-y)dy 



So we must find an upper bound for r n (j — k). Let us note that the function ">p(h) = 
/i _7 / 2 L(/i) is regularly varying of index —7/2 < 0. Therefore, the Karamata theorem of 
uniform convergence (see Bingham et al. [7], Theorem 1.5.2) yields 

ip(zh) _ 1 



lim 

h^o+ ip(h) z^l 2 

uniformly in z G [3, +oo[. 

As a consequence, there exists K > (which depends only on L and 7) such that for 
n large enough, 



Vz > 3, tp [ - ) < Kip 



(33) 



However, we have 

r n (j-k) = n^- 2 
Therefore, (17) yields 



j+1 



du / dv 

'u-l 



k+l 



dx 



r„(J -k)< Kn^-'L ( - ) p l/2 (j - k) < K 



Consequently, wc have 



djk -C[ -,- )r n (j - k) 
n n 



< Ke Tl 



x ip{(v-y)/n) 
-1 {v-yy+ti/ 2 ) 

n<- 2 L{l/n) 
(J -k- 2) 2 +(t/ 2 ) ' 

n>- 2 L{l/n) 
(j-k- 2) 2 +(^ 2 ) ' 



dy. 
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and because C is bounded, 



C( J -,-) r n (j-kf 



<Ke n 



n 2 T- 4 L(l/n) 2 
{j-k- 2)4+7' 



Using the same perturbation argument as in Begyn ([5], pages 10-11), we can check that 
it is still true whenever j = 1 or k = 1. Consequently, 



,3-27 



L(l/n)= 



-2n-l-f 



i — 1 — Z 



E E <W-E^(0 2 E c 



1=3 k=l 
n — 2 n—l—l 



1=3 



k=l 



k + l k 



n n 



< 



A " e «-E E ( /_ 2) 4 +7 <^Et7^ 



(l- 2) 4 +7 



;=3 fc=i v ' 1=3 
and the right-hand side is convergent because 4 + 7 > 1 . This yields 



lim 



,3-2 7 



n->+oo L(l/n) 2 
Moreover, 



<n— 2n— 1— i n—2 

E E <W-E^(0 2 E c 

U=3 fe=l Z=3 fe=l 



(34) 



lim I y C(— ,-) = lim - Y C{-,-) = I C(x,x) 2 dx 

— >4-r>n 71 ' J \ 71 71 I n. — >-t-oo 71 f 1 \ m n I /„ 



l-l-l 



k=l 



k + l k 



l-l-l 



ii hoc /( 



k=l 



k k 



n n 



thanks to the uniform continuity of C. 

In addition, because L is slowly varying, the theorem of dominated convergence yields 
(using inequality (33)) 

r„(0 



lim 

n->+oo rCi- 2 L(\jn) 



p~t(l). 



Hence for all I > 3, 



lim 



rw+oo rt 2 7— 4 L( 

and (using (33)) we can check that VI > 3, 



.(<) 



l/n) z n z — ^ \ n n J " 



l27-4£,(l/ n ) 



A" 



- 2~l 4 +7' 



(i-2) 



Consequently, the theorem of dominated convergence for series yields 



,3-27 



i-i-/ 



i^+oo L(l/n)'' 



E'«(o 2 E c 



k=l 



k + l k 



+ DC 



C(x,x) 2 dx^p 7 (/) 5 



i=3 



Quadratic variations of Gaussian processes 
With (34), we obtain 
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n-2 n-l—l 



J&o 1(17^ E <4 - „ i™. r^/nf E E 



j-fc>3 



.Step 2.2. Let us prove that 

„3-2 7 

lim 



n — >+oo 



L{l/nf 



«1 +oo 

/ C(s, a:) 2 da: ^^(0 
Ja 1=3 



J2 d% = Pl (2) 2 / C(x,x) 2 dx 



l<k<j<n-l 
j-k=2 



With the perturbation argument of Begyn [5] , we can check that 



i k ~\~ 2 k i 

Mfc+2,fc - CI , - r„(2) 

n n 



<tfe„^- 2 L - p 7/2 (2) 



and, consequently, 



dl+2,k C 



fc + 2 fc 



>(2) : 



<Ke n n 2 ^L - p 7/2 (2) 



Therefore, using the same arguments as in Step 2.1, we have 



lim 

n^+Ki L(l/n) 2 



E = 7777^2 E 



l<fc<j<ra— 1 
J-fe=2 



n— >+oc L(l/n) 2 ^— ' 

r„(2) 2 1 



lim 



i — 3 /fc + 2 fc 



E^ 



ra— >+oo 7i 2 T 4 L(l/n)n^— ' \ n ' n 
p 7 (2) 2 / C(x,x) 2 dx. 



5"tep Let us prove that 



n 3-2 7 _ /•! 

lim . , - 2^ d j7 . = / ff (x) 2 dx. 



n->.+oo L(l/n) 5 



l<fe<j'<n-l 



n 



This is a consequence of the assumption (23). 
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Step 2-4- Let us prove that 



n 3-2 7 rt 
lir ? 7777^2 E d lk = / 5o(a;) 2 da;. 



fc=l 



(37) 



This is a consequence of 



lim sup 

h^0+ h<t<l-h 



L{h)h 2 -~< 



9o(t) 



:(). 



which comes from the assumption (22). 

The preceding four steps imply (31). Let us remark that (30) and (31) yield that 



Step 3. To prove Theorem 2, we use the decomposition 

T n - \fnb n = T„ + ET„ - \fnb n . 



Let us prove that 
We have 

E(T„) - \fnb n = \fn j 
On the one hand, 



lim \ET n - y/nb n \ = 0. 

n — ►+oo 



L(l/n) 



E dfcfc ~E/ gi(x)dx(f>l - J 



l/n) n f—f V n / z — ' V n 



L(l/n) 



fe=i 



i^i n — 1 



fe=l 



< v 7 "- sup 

fc=l n-1 



Hfe V - /AS 

Lfl/n) ^ ff Hn 



n7- 2 L(l/n) 



Therefore, the assumption (22) yields 



lim \fn 



n — >+oc 



T(l/n) 



fe=l i=0 

Moreover, if we choose < i < g— 1, we have 



:-r^ 4fe - - v f -\ y^gJ- 

' ; fe=l i=0 v 7 fe=l v 



1 fl 



n \n 



Eft 

fc=i 



Pi (a;) dx 



< K(j) 



iy i 



(38) 

(39) 
(40) 



(41) 
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due to the fact that is bounded and Lipschitz on ]0, 1[ (see the Proof of Theorem 1). 
Furthermore, 



T) \ T) I — ' \ n I \ T) 



fe=l 



g q (x) dx 



< KM - 



knowing that g q is bounded and 1/2 + aq-Holderian on ]0, 1[. Consequently, 



lim v» 

n — >+oo 



■El 

i=0 



Yj9i 

fc=l 



gi(x)dxcj) 



(42) 



To finish the proof, just note that (40) is a consequence of (41) and (42). 

Next, by combining the Prokhorov theorem, all the preceding steps and the Slutzky 
lemma with (39) and (40), we get (24). □ 

In the sequel, we consider estimators of some functions of the parameters 7, which 
are constructed with both second-order quadratic variations V n (X) and V2 n (X). So the 
preceding theorem is not sufficient to prove the asymptotic normality of the estimators. 



2.3. Bivariate central limit theorem 

The next theorem will be useful to prove the asymptotic normality of our estimators. 
We define the following constants, which appear in the asymptotic covariance of the two 
quadratic variations: 

a? =2/ g(t) 2 dt + 4p 7 (2) f g (t)C(t,t)dt + 4p 7 (3) [ g(t)C{t,t)dt 



pi +00 
+ 4 / C(M) 2 di5> 7 (Z)p 7 G-2), 

Jo 1=4 



4cov = 4/ go(t)g(t)dt + 4p 7 (2) / g(t)C(t,t)dt 



pi +°° 

+4/ c(t,t) 2 dtj2p 7 (i) Pl (i-i 
J ° 1=3 



and 



a 2 = 3o- 2 + a 2 + 4cr| 



(43) 
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Theorem 3. We use the same assumptions as in Theorem 2. Then we have 



%^V 2n (X) So g (x) dx - ELi Jo 9i{x) dx ■ 
where the 2x2 matrix £ is defined by 

£ 



JV(0,£), 



a 2 2^-V 2 

-2„2 „2 



2T- 2 ^ a 2 /2 



Proof. We set for A, fj, G R, 



We begin by showing that when A, > 0, 

VnlSnfAj/jJ-AV / ffi(x)da;0( - J / ffi(x)da; 



First we must prove that 



We have 



lim nVarS^A,^) = AV +fJ?— + 2'>- 1 Xna 2 . 

n — >-j-oo 2 



„2-2 7 (2n) 2 ~ 2 '> 
VaxS n (\,n) = A 2 — — Var F„ (X) + M 2 /. . .„ Var V 2n (X) 



L{l/n) 2 ^ " " L{l/{2n)f 

2 1 -7 r ,2-27 

+ 2X » tim U7T775 ^ Gov(y„(X), F 2n (X)). 
L(l/n)L(l/(2n)) 

In the Proof of Theorem 2 we showed (31): 

lim F7i / v> VMV n {X)=a 2 and lim / / V8xV 2n {X) = (j 2 . 
n^+oo L(l/n) 2 n^+oo L(l/(2n)) 2 

We must compute the term Cov(V„(X), V 2n (X)). We have 
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Therefore, 



V n (X) = Y J \^xg%? + {Axg n \f + A(Axg n) ) 2 
fe=i 

+ 2AX^ 1 Axg n \ + 4AXg%AX™ + iAX^AXg n \ 

Moreover, 

2n-l 

V 2n (X)= ^(Axf n) ) 2 . 

To simplify, we will set the notation AXi for AX^ 2n ^ and di p for d?" . To compute 
E[V r „(X)V2n(^)], we use the Isserlis formulas (see Isserlis [14]), which yield 

Cov[(AX/) 2 , (A^-) 2 ] = 2(E[AXiAXj}) 2 = 2d% 

and 

Cov[AXiAX pi {AX.jf\ = 2E[AX i AX i ]E[AX p AX J -] = 2d lj d pj . 
So we can check that 

6 



Cov(V n (X),V 2n (X))=J2s t , 



i=l 



with 



ra-12n— 1 n— 12n-l 

fe=i j=i fc=i j=i 

n-12n-l n-12n-l 

k=l j = l k=l j = l 

n-\1n-\ n-\1n-\ 

<S , 5 = 8^^ d2k+i,jd2k,j, 5"6 = 8 d-2k-i,jd2k,j- 

k=l j=l k=l j=l 

However, using the same techniques as in the Proof of Theorem 2, we obtain the formulas 

2 1 — 'n 3-27 
lim \ i~ ? -t ,, n \ \ Si = 2 7 ~ 3 cr 2 

L(l/n)L(l/(2n)) 

2 1 — r n 3 ~ 27 
L(l/n)L(l/(2n)) S4 = T '~ a ^ 



2 1 — 'n 3-27 

lim ; — ; ttSr = 2 7_1 (T 2 

L(l/n)L{l/(2n)) 
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and, likewise, 



2 l-7 n 3-2 7 

lim 

n^+oo L(l/n)L(l/(2n)) 

2 i -7 n 3-2 7 

n^+oo L(l/n)L(l/(2n)) 

2 l-7 n 3-2 7 

™-+oo L(l/n)i(l/(2n)) 



5 2 =2T-V, 
S 3 = 2T-V, 

^6 = 2 7 1 cr| . 



Hence, 



lim 2 1 -Tn 3 " 2 T 
ni+oo L(l/n)L(l/(2n)) 



Cov(y n (X),y 2 „pO) = 2^- 2 ^ 



(49) 



Therefore, (47) is a consequence of (48) and (49). 

Now we apply the Lindcberg central limit theorem to S n (X,fi) in the same manner as 
in Theorem 2. We set 

S n (X,fi) =S n (X,fi) -ES n (X,n). 

Because A > and \i > 0, we can consider ^(A, \x) as the Euclidean norm of the Gaussian 
vector (d; 1 < i < 3n — 2) : 



l<i<n-l, 



L(l/(2n)) 
Therefore, Cochran's theorem yields 



5„(A jM ) = X;r J> (y«) 2 , 



with a„, Tj n , r* and 5"^ as in the Proof of Theorem 2. This yields 



5„(A > /i)=£r iin [(F n W) 2 -l]. 



Also, we notice that 



3n-2 



t* < max V |E(ZiZ ? -)|. 

1<7<3ti-2 ' J 



i=l 
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Moreover, if 1 < i < rt — 1 and 1 < j < 2n— 1, 



j(2") , j(2n) „ ,(2n) 

u 2i+l,j u 2i-l J ' zu 2i,j ' 



n-1 2n-l 



7"n - K-fj—j-r ( A max V|4 n) |+M max V \d { ^ n) \ 

L(l/n)\ l<j<n-l^ J l<j<3n-l J 



i=l i=l 
2n-l 



(2n), 
»,2j I 



l<7<n— 1 J J 

i—1 

n-1 

V max E(l4-?i J + l4*-U + 2|4-Jl) 

l<7<2n— 1 * — ' ,J ,J >J 



n-1 2n-l 

<K^j—\ max V|4 l) |+ max V |dg n) ' 

L(l/n) \ l<j<n-l ^— ' iJ l<j<2n-l ^ 3 

' \ z— 1 i—1 

Therefore, with the arguments of Begyn [5], we can checked that 

r* < KrC~ 2 L 

Thus 



t* K 
~ v/VarS„(A, M ) v 7 ^ 



Using the Lindcbcrg central limit theorem, we obtain 

(A, A1 )-E5„(A,/i) (£) 



• v /Var5 n (A,/x) 



AA(0,1). 



Hence, as in the Proof of Theorem 2, we establish the convergence announced in (46). 
With a generalization of the Cramer-Wold arguments, based on the properties of the 
Laplace transform that can be found in Istas and Lang ([15], page 431), we ascertain 
that the Laplace transform of the vector 



gKW - Si 90 (t) dt - E?=i Jo 9i(t) dt ■ 0(1) 



( ^rV 2n {X) - £ 50 (t) dt - ELi Jo M*) d * ' <!>& 



tends on R 2 to the Laplace transform of a Gaussian centered law with covariance matrix 
E. This proves the result. □ 
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As we stated in the Introduction, we want to apply these results to estimate the 
parameters of some fractional processes. To explain how to use these results, we apply 
them to the FBM, even though the consequences are not new in this case. We refer to 
Sections 3 and 4 for original consequences. 



2.4. Application to fractional Brownian motion 

We study the example of the FBM B H . We can check that the theorems of Section 2 
can be applied with 7 = 2- 2H, L(h) = 1, q = 0, g (t) = 4 - 2 2H and = h. 
By applying Theorem 1, it follows that almost surely (see Cohen et al. [10]) 

lira n 2H - 1 V n (B H ) = 4-2 2H . (50) 

n — >-\-oo 

Next, if we apply Theorem 2, then we obtain (see Coeurjolly [9]) 

^{n 2H - l V n (B H ) - (4- 2 2ff )) ^M(0,4bmm) (51) 

with 

°Ibm,h = 2(4 - 2 2H ) 2 + (2 2H+2 - 7 - 3 2H ) 2 

+ (2H) 2 (2H - l) 2 (2H - 2) 2 (2H - 3) 2 \\p 2 -2H || 2 , (52) 

because, when computed, C(s, t) = -H(2H - 1)(2H - 2){2H - 3) and g(t) = (2 2H + 2 - 
7-6 2H )/2. 

Because of Theorem 3, we get 

^V n {B H ) (4 - 2 2H ) \W Kf(() y , (w 

\2H-1 V (y\ (a 2H\ ) JV (0, i^FBM.H ) (53) 



(2n) 2H - 1 V 2n (X)-(4-2 2H ) 



with 



J FBM,if 



2 9- 2H rr 2 
"FBM,H z a *,FBM,H 

z °*, FBM, if z °FBM, if 



(54) 



and 

^?,cov,FBM,H = (2H) 2 {2H - l) 2 (2H - 2) 2 {2H - 3) 2 J2p2-2h(1) P 2-2h(1 - 2) 



1=2 



+ ±(2 2H+2 -7-3 2H ) 2 , (55) 

+ 00 

4cov,fbm,h = (2H) 2 {2H - l) 2 (2H - 2) 2 {2H - 3) 2 ]T p 2 . 2H (l)p 2 - 2H {l - 1) 



1=2 

-,2H\/r,2H+2 n n2H\ 



+ 2(4-2^)(2^+"-7-3^), (56) 

CT *,FBM,if = 3o"FBM.H + ^l.cov.FBM.ff + 4c>2,cov, FBM, if • (57) 
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The 5 method yields that the statistic 

1 , V2n(B H ) . .J 

^ = 2- log K(^y (210g2) 

is a strongly consistent estimator of H and that (see Coeurjolly [9]) 

/ 0—2 r\2 — 2H 2 \ 

- H) _» A^O, 4(4 , 22g)lQg2 J • 

3. Two- parameter fractional Brownian motion 

Two-parameter fractional Brownian motion was introduced by Houdre and Villa [13] 
as an example of a quasi-hclix. Two-parameter fractional Brownian motion B H,K = 
{B^' K ;t £ K} is defined, for H e]0, 1[ and K e]0, 1], as the unique continuous centered 
Gaussian process with covariance function 

\/s,t £ R i? H ^( s ,t) = Cov(B?> K ,B?> K ) = ±((s 2H + t 2H ) K -\s- t\ 2HK ). 

The process B H,K is i?X-self-similar, it has a critical Holder exponent equal to HK in 
the sense of Adler [1], it is, for K = 1, the standard fractional Brownian motion and it 
has stationary increments if and only if K = 1 . 

We refer to Houdre and Villa [13] for the proofs. In their paper, they introduced the 
process B H ' K under the name bifractional Brownian motion. We suggest here to call 
it two-parameter fractional Brownian motion, because it is a monofractal process (the 
pointwise Holder exponent of its trajectories is a.s. constant) and the term 'bifractional' 
may imply that it is a multifractal process with two values for its pointwise Holder 
exponent. 

Theorems 1 and 2 yield estimators of H and K, the relevant quantities in the study of 
B H ' K . To determine the estimators, that we assume that we dispose of the observation 
of one path of B H ' K on the interval [T l5 T 2 ], where T U T 2 £ E, T x < T 2 . This process is 
considered to be indexed by [0, 1]. Therefore, we introduce the process Y defined by 

Vie [0,1] Y t = B H - K { T {t)) withT(t) = (T 2 -T 1 )t + T 1 . 

We obtain a new process Y, which is centered, Gaussian and has covariance function 
r H - K {s,t) = Cov(Ks, Y t ) = R H ' K (r(s),r(f)), and we dispose of the observation of one 
path of Y on [0,1]. The results of Section 2 can be applied to Y under the condition 
[T X ,T 2 ] C]0,+oo[. 

3.1. The results 

First we study the almost sure convergence of the second-order quadratic variations. 
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Proposition 4. We have when n — ► +00 

4 



* HK - l V n {Y) ^ 4 J-i {T2-T,f HK . (58) 



Next we study the weak convergence. 
Proposition 5. We have, when n — > +00. 

V^(n 2 ^-V„(F) - ^^(T 2 - T X ) 2 ^ 

where <Tp BM was defined in (52). 

As in the case of the FBM, we can deduce an estimator of HK. 
Proposition 6. The statistic 

^„ = 5-log^^(21og2)- 1 (60) 

is a strongly consistent estimator of HK and when +00 , 

^(ffA ri -FA)^A^0, 2 *+ 1 (4-2 2 ^)log2 (T2-T1) J, (61) 

where OpBM //if (T * fbm irif were defined in (52) and (57). 

The quantity 7JA" is relevant in the study of B H ' K , but it does not characterize the 
law of this process. For this characterization, we need to know both parameters H and 
A. A refinement of the previous results enables us to construct strongly consistent and 
asymptotically normal estimators of these quantities. 

Proposition 7. The statistic 

1 / n 2HK n -l \ 

K n = 1 - — log ^ ^- V n (Y) (62) 

log 2 \(4-2 2HK »)(T 2 -T 1 ) 2HK n / 

is a strongly consistent estimator of K and when n—t +oo, 

^ ), (63, 
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where Cp BM HK was defined in (52). Moreover the statistic 

S„ = Si ,64, 
A',, 

is a strongly consistent estimator of H and when n — > +oo , 

Vn-(H n -H) HmU &0??l rA (65) 



with 



and 



2 H 2 1 2H 



CT FBM, HK 

(4 - 2 2HK ) 2 log 2 2(T 2 - Ti) 4HK ' 



o 2 _ 9 2~2HK 2 

V2 — z 



4(4-2 2ff *")log2 

0-2HK 2 2 

z a *,FBM,_f/X "FBM.ffJi' 



2(4 - 2 2HK ) 2 log 2 2(T 2 - T X ) 2HK ' 
where a 2 FBM was defined in (57). 

3.2. Proofs of the results for two-parameter FBM 

Proof of Proposition 4. We apply Theorem 1 to the process Y. We need to show only 
that assumptions 2(b) and 2(c) (in Theorem 1) are satisfied (the other assumption is 
obvious). 

For assumption 2(b), it is clear that the derivative g s 2 g t 2 (s, t) exists on ]0, l] 2 \ {s = t}. 
Moreover, we can check that, Vs, t € ]0, l] 2 \ {s = t}, 

d 4 r H K 2HK(2HK - 1)(2HK - 2)(2HK - 3) 2ra , 2HK _ 4 

d^W^ 2^ (T2 ~ Tl) |S " M 

+ (T 2 -T 1 ) 4 ^(r(s) 1 T(t)) 1 (66) 

where i/j(t(s), r(t)) is continuous on [0, l] 2 . Therefore, the assumption 2(b) (in Theorem 1) 
is satisfied with L(h) = 1 and 7 = 2 — 2HK . 

For assumption 2(c) (in Theorem 1), computations yield 

(ffo^V^XM) 4-2 2 ^ T , 2HK .e t {h) 
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and we can check that £{(0) = e' t (0) = e'/(0) = e\ (0) = 0. So that Taylor formula yields 

\/h<t<l-h e t (h)= [ h ( - k ^. X)3 ef\x)dx. 



Therefore, we have 



.3! 



sup sup \e y t'{x)\ = 0(1), 

h<t<l-hO<x<h 



which yields 



sup 

h<t<l-h 



(6>?o5Zr H ' K )(t,t) fm m , 2HK A-2 



2HK 



(T 2 — Ti) z 



0{h^ 2HK ). (68) 



h 2HK \-* 2 K ~ X 

Therefore, the assumption 2(c) (in Theorem 1) is fulfilled with 

ffo(i) = i ^zi-(T 2 -T 1 ) 2 ^. (69) 
Consequently, we can appply Theorem 1 to Y and obtain (58). □ 

Proof of Proposition 5. We apply Theorem 2 to the process Y. As in the Proof of 
Proposition 4, we need to show only that the assumptions 2 and 3 (in Theorem 2) are 
satisfied. 

For assumption 2 (in Theorem 2) the previous proof showed formula (66), which yields 
that for all s, t e]0, l] 2 \ {s = <}, 

d 4 r H - K 2HK(2HK - 1){2HK - 2)(2HK - 3) vihk i ,^hk-a 

+ (T 2 -r 1 )V(T(s),r(t)) J 

where ip(T(s),T(t)) is continuous on ]0, l] 2 . Therefore, the assumption (21) in Theorem 
2 is satisfied with L(h) = 1, 7 = 2 - 2HK and 

r( ,, 2HK(2HK - 1)(2HK - 2){2HK - 3) 2HK 

+ (T 2 - 7i) 4 | S - i| 4 - 2 ^V(r( S ), r(t)). (70) 

For assumption 3 (in Theorem 2), formula (68) of the previous proof shows that the 
assumption 3(d) (in Theorem 2) is fulfilled with q = 0, go(t) = 8 ~ 2 2 k an d «o = 1/2. 
Moreover, we can check that 

(ff o gfrgfXt, t + h) 2 2g *+ 2 - 3 2 ^ - 7 _ T s 2 i?if , Vt0})_ 

foZHK 2 K h 2HK ' 
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With the same arguments as those used for et(h) in the previous proof, we obtain 

Sup \T] t (h)\ h ^ + G(h 4 ). 
h<t<l-h 

This shows that the assumption 3(e) (in Theorem 2) is satisfied with 

2HK+2 _ o2HK _ 7 

= w (Ta Tl)2HK - (71) 

Consequently, we can appply Theorem 2 to Y and obtain (59). □ 
Proof of Proposition 6. We apply the 5 method with the C 1 function 

)= i_ iog(y/s) 

nX,V) 2 21og2 

to the convergence announced in (44) to yield the result. □ 

Proof of Proposition 7. First we establish a refinement of Proposition 4. Because of 
(68), we have, for all a G ]0, 1[, 



sup 

h<t<l-h 



(#oJjr*-*)(M) (T T ,2HK%-* HK+l 

{J- 2 - ilj 



Jj2HK ^ ->-J 2 K 



h =° + o(fr Q ). 



Therefore, the assumption 2(c) of Theorem 1 is fulfilled with g = 1, # (i) = {T 2 -Ti) 2HK : 
4 ~^_! , <7i(f) = 0, 0(/i) = /i and i^i = a. It yields that almost surely 

A _ 2 2HK / 1 



i V n (Y) = (T 2 -Ti) x +o 



Taylor expansions yield that almost surely 



and 

n 2HK n -i n _ d _ oo ^logn 

n 2HK-\ ~ t + O 

With a = 3/4, we obtain that almost surely 
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In addition, we have (44): 

-,2HK-1m r v \ 1-2'"^ (m m\2HK \ , „, / (m \4HK 



If we apply the 5 method with the C 1 function 



x 



f{x,y) = I 1 _ logfa/s) , 

\2 21og2 / 

then the Slutsky lemma and (72) yield that there exists a 2 x 2 real matrix A such that 



TO-^ (£) kr ( n (T 2 -T 2 ) 



\iHK 



v/i| ( 4-2^K„)( T2 _j. 1 )2H K „ 2 >AM0, n2(Jf — 1) 

By again applying the (5 method with the C 1 function 

/ -, _ logic 

we obtain 

Kn-K \ ^, A f(Q} T ^- T ^ HK , 



with 



HK n -HK 1 2 2 ( A '-! 



??3 m 



This proves (63). A final application of the S method with the C 1 function 

f{x,y) = - 
x 

yields (65). □ 



4. Anisotropic fractional Brownian motion 

Let c(sN*. Let H : M. d — > ]0, 1[ be a Borelian function that is homogeneous of degree zero, 

V£eR d ,VAeM\{0} H(X£_) = H{£_), 

that can be identified with an even function from the sphere S^ 1 into K that we denote 
H as well. We assume, moreover, that H takes its values inside the interval [H_, H] C ]0, 1 [, 
with H_ = ess inf H and H = ess sup H. 
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We define the anisotropic fractional Brownian motion (AFBM) with directional Hurst 
index H, denoted A^ H \ by the harmonizable representation formula 

VueR d A^(u) = J Rd ^§^dW(a (73) 

where (•,•) is the canonical scalar product and W is a complex random measure in the 
sense of Samorodnitsky and Taqqu ([16], 325-328). It is a Gaussian field with stationary 
increments. Bonami and Estrade [8] showed that A^ H ^ has a critical Holder exponent 
equal to H_ . Moreover, they showed that the field A^ is locally asymtotically self-similar 
(l.a.s.s.) of order H at any point of M. d (see Definition 8) if and only if Leb({H(9) = 
H}) > 0, where Leb(-) denotes indifferently the Lebesgue measure on M. d or the Lebcsgue 
measure on 

Let us recall the definition of the l.a.s.s. property (see Benassi, Jaffard and Roux [4]): 

Definition 8. Let (3 > 0. A process {X u ;u g M. d } is locally asymptotically self-similar 
(l.a.s.s.) of order (i at point Uq G R d if the finite- dimensional distributions of the process 

X(u + Xu) -X(uo) 
\P 

converge to the finite- dimensional distributions of a non-zero Gaussian process when 
A — > + . The limit process is called the tangent process at point uq . 

Our purpose is to identify the function H when we consider one realization of the field 
A( H ' . For that we apply the theorems shown in Section 2 and restrict the field to some 
segment of R d . 

To simplify the computations we assume next that d = 2. Note that in this case we can 
identify H with an even 7r-periodic function on R. We consider one realization of A^ H > , 
which is observed in axes denoted by Oxy. We assume too that these axes of observation 
are equal to the axes of definition of A^ . 

Let [A, B] be the radial segment of length L <G ]0, +oo[ such that the distance between 
O and A is equal to Le (with e > 0) and the angle between [^4, B] and the axes Ox is 
equal to uj £ [0, 2n[. See Figure 1 for more details on the geometry of the problem (and 
note that the angle are oriented anticlockwise). 

We use the following parametrization of the point u = (m, U2) E [A, B] : 

u\ = L(t + e) cosw, 

, x • (74) 

112 = L(t + e) sintj, 

where t goes over the interval [0,1]. Next we consider the restriction of the field A^ H ' 
to the segment [A, B] with the parametrization (74). Hence we obtain a new process Z 
indexed by t £ [0, 1]: 

(l0 /L(t + e )co8 W 
1 \L(t + e)smw 




It is clear that Z is a centered Gaussian process with stationary increments. 

To apply the theorems of Section 2, we must compute the covariance function R of Z 
that is given by the following lemma. 

Lemma 9. The covariance function of the process (Z t )t^[o 1] is given by, for all s,t<E 
[0,1], 

R(s,t)=4 f A(6)[\ s + e\ 2H ^ + \t + e\ 2H ^-\ S -t\ 2H W]d6, (75) 
Jo 

with 

We [0,2*] A(9) = gg^p ^W|cos(g-a;)| 2g W, (76) 
where for all d €E N* and for all H G ]0, 1[, 



/ 71^+1)^(^+1/2) \ 1/2 



C(d,H) - [ HTi2H)siniH7r)riH + d/2) ) (77) 
cme? r denotes the Euler gamma function. 

Proof of Lemma 9. First we compute the variogram v of the field . For all u£R 2 , 
we have 

v(u) = ±E[(AW(u)) 2 } = l I C(l,H(y)) 2 \(u,y)\ 2H ^dy 



i r^rc 



s 1 

C(l, £T(0))|m cos(0) + u 2 sin((9)| 2H(e) d0. 
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Because the field A^ has stationary increments and is vanishing a.s. at the origin, its 
covariance function is given by 

Vu, u'el 2 , Cov(AW (u), A (h > (u')) = 2(v(u) + v(u') - v(u - u')). 

If we take u, u' in the segment [A, B] with u parametrized by s and u' by t, we obtain 

R(s,t) = Cov{AW(u),AW(u')) 

= 2 r k{6)[\s + e\ 2H ^ + \t + e\ 2H ^ -\s- t\ 2H ^} dO. 



o 



Because the functions A and H are 7t-periodic, we obtain (75). □ 

By applying theorems of Section 2 to the process Z, we are able to estimate H_ , the 
Holder critical exponent of . We distinguish two cases. 



4.1. The l.a.s.s. case 

In this subsection we assume that 

Leb{H(6) = essinf H} > 



and we use the notation 



vi6]o,i[ Jh = 8 [ A(e)i {H{e)=m de. 

Jo 



(78) 



Proposition 10. We have, when +oo, 

n m-i Vn ( Z )^(4-2 2 Z)J K . (79) 

For the central limit theorem we must study the asymptotic of 

h -> 8 [ n A(6)(4 - 2 2H W)\h\ 2 W e ^h {K<H(e)< H +1/4} d9 
Jo 

when h — > + . For this purpose, we need to sharpen the assumption on the function H. 
We get two kinds of central limit theorem. 
Case I. Assume that 

Leb({ff < H{9) <H+l/4}) = 0. 



742 A. Begyn 

Proposition 11. We have, when n— > +oo, 

V^(n 2E -~ l V n {Z) - (4 - 2 2 ^) Jh) H Af(Q, JWfbm.h)- (80) 



From these results, we can deduce a strongly consistent estimator of H_ that is asymp- 
totically normal. 



Corollary 12. The statistic 



is a strongly consistent estimator of H_ and 



o 2 _ 2 2-2H 2 



Vn(H n -H) >Af[ 0, Jh 4( 4 _ 2 2g )log2 ) . ( 82 ) 



where cr| BM H and a\ FBM H are as defined in (52) and (57). 
Case II. Assume that 

Leb({F < H{6) < H + 1/4}) > 0. 
Proposition 13. We have almost surely 



n^^V n (Z)'^(4-2^)J K +^-j +°{ (l) {n)) (83) 
with 

m = 8 f A(0)(4 - 2^)\h\^-^h {K<H{e)<K+1/i} AO. 
Jo 

Moreover, when n — * +00, 



>[ljj^Af(0,jy 2 FB M >M )- (84) 

As in Case I, we can deduce a strongly consistent estimator of H_ that is asymptotically 
normal. 



Corollary 14. The statistic 



^l-^WT (2loe2rK (85) 
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is a strongly consistent estimator of H_ and 

r frr ffx] A4-2 2 g)JH + ^(l/(2n)) y oi 

/ _ 02-2/^2 x 

( £ ) irl n t FBM,iJ z — U *,FBM,H \ 

— 4^—2^)1^2 => (86) 

Remark. In general, we have 

(4-2^)^ + 0(1/^ J (21 ° g2) =°' 

but the limit 

^ (4-2^)Jg + 0(l/(2n)) y oi 
lmi Vnlog — — n ,„- 77—7 — r (2 log 2) 

does not exist in general. 

Exam/pie. We assume that 

H(6)=Hl ]QM (e) + Hl ]eoM (9) 

with H_<Tl and < 6» < 7r. 
If If > if + 1/4, this is Case I. 

If H = H + 1/4, this is Case II. We can check that = (4 - 2 2 ^ +1 / 2 ) Jff+i/ 4 /i 1/2 , 
which implies 

^ (4-2^)J g + ^(l/(2n)) ^ oi (4-2 2 ^+V2)J 1 ^2 

hm Vnlog — — — — — — — — (2 log 2) L = 



«-+oo v b V (4-2 2 ii)JH + 0(l/n) / v b ; (4 - 2 2 ^) J K 2^2 log 2' 

Therefore, the Slutsky lemma yields that \/N n (H r — H_) converges in law to a Gaussian 
random variable with mean equal to the opposite of the right-hand side of the previous 
formula. 

HH<H + l/4, then this is also Case II and we have <f>(h) = (4 - 2 2 ~ S )J 7T h 2( - 7I -S-\ 
which implies 

Jffoo^H (4-2 2 il)JH + 0(l/n) J (2l0g2) ^ + °°- 



In this case, the bias term becomes infinite, so in practice it cannot be ignored. 
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4.2. The non-1. a.s.s. case 

In this subsection we assume that 

Leb({H(9) =essinfF}) = 0. 

Next, we add the following assumptions: 

1. H is of class C 2 on R and there exists a unique point 8* £ [0,7t[ such that H{0*) = 
H = w£H. 

2. H is decreasing on ]£>* — 7t/2, #»[ and if is increasing on ]6>*,6>* +7t/2[. Moreover, 
#"(0*)>O. 

3. The angle w is such that to ^ 6* + n/2 (mod7t). 

These assumptions are classical when the Laplace method is applied. It is well known 
that they can be weakened, but we skip the technical details (see Dicudonnc [11], page 
125). 

Proposition 15. For all to £ R 2 , the finite- dimensional distributions of 



( _ log£) v,^±4z^M ;i 



eR^| (87) 

converge when e — > + to the finite- dimensional distributions of the process 
C{1,H)( 1 \ 1/4 Dff ,^ ^.^ TO2 " 



BH Pe ,mt£ 



where p$ t (t) is the orthogonal projection of t on the straight line {rc l8 ' , r £ R} and B— 
is a standard FBM of Hurst index H . 

Consequently, the field A^ H ' is no longer 1. a.s.s. and we have shown that a normalization 
with a logarithm factor yields a non-trivial limit field. We will see in the sequel that we 
have the same behavior for the second-order quadratic variations. 

For this section we use the notation 



H"{6*) 

Proposition 16. We have, when +oo, 



^- 1 ^irlV n (Z)^(A-2 2 ^)Gg,. (89) 



Now we want to show a central limit theorem. For that we must strengthen assumptions 
on H that were demonstrated in the beginning of Section 4.2: 
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1. H is of class C°° and there exists a unique point G [0,7r[ such that H{9*) = H = 
mfH. 

2. H is decreasing on ]0* — 7t/2,0*[ and 7? is increasing on J^*,^* +7t/2[. Moreover, 
#"(0*)>O. 

3. H is analytic in the neighborhood of the point 

+ OC 



H(6) = H + J2 H i 



i=2 

4. The angle u> is such that uj ^ 6* + n/2 (mod7t). 

The assumptions imply that the function A is also analytic in the neighborhood of the 
point (9*: 

i=0 

In this case, we use the extended Laplace method (see Wong [17]). 
Proposition 17. For all ggN, q > 1, we have, almost surely, 

n- L sf\ognV n {Z) = 16> — -j- — t + o[- — ^ , 

(logn)v^ y^logn)?/^/ 

where the coefficients 0{ can be expressed in terms of Hi and Ai . 

Remark. We can give that explicit forms for the coefficients a.i (see Wong [17]). For 
instance, the first two coefficients are given by 

4-2 2 -^ 

a = __G e ,, (90) 
16V^ 

2 / A'(fl.)(4-2 2 £) H^(6*)A(6*) \ 
1 H"(6*)\ 2 3i/"(^) / 1 ' 

Proposition 18. We have, when +oo. 

/ rn f , x 2(H(9)-H) x 

V^loi^^-" 1 ^^)-^ A(0)(4-2 2H W)^-J l {ffW < S+ i/4}dflJ 

^^(0,G^<r§ BMii£ ). (92) 

Theorem 3 and the Slutsky lemma yield that when n — > +oo, 

/— — / n^- 1 K(^-/ o 7r A(0)(4-2 2 ^))(i) 2 (^)-^l {H(e)<ii+1/4} d0 \ 
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(C) 



AA(0,G^S fbm ,h). (93) 

As in the previous case, we obtain an estimator of H_ that is strongly consistent and 
asymptotically normal. 



Corollary 19. The statistic 



^I-^Wj^" (94) 



is a strongly consistent estimator of H_ and 



-^^{^ 4(4 -2^) log 2 ' 



with 

Hh)=8^]^h f m^-^ H(e) )\h\ 2[H(e) - K) l { H { e)<H + i,A } ^ 
Jo 



Let us note that because Z is Gaussian and has stationary increments, we can apply the 
results of Istas and Lang [15] to estimate H_. To do so we must assume that Leb({H(9) = 
H}) > and H_ > 3/4. Moreover, we need observations of Z along an infinite interval, 
which is not the case in our assumptions (t £ [0, 1]). In this sense, we have improved the 
result of Istas and Lang [15] in the case of the AFBM. 

Let us note that the estimation of the function H was performed by Ayache et al. [2] 
and Bierme [6]. 



4.3. Proof of the results in the l.a.s.s. case 



Proof of Proposition 10. We must check that the assumptions of Theorem 1 are 
satisfied. For assumption 2(a) (in Theorem 1), note that the functions H and A arc 
bounded functions that can be deduced by the Lcbesgue theorem of continuity under the 
symbol integral. 

For assumption 2(b) (in Theorem 1), we must compute the derivative q§2-§£i- For 
the same reasons as above, the Lebesgue theorem of differentiability under the symbol 
integral shows the existence of this derivative on ]0, l] 2 \ {s = t} and yields the formula 
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A(0)2#(0)(2#(0) - 1)(2H(6) - 2){2H{6) - 3)|s - t\ 2{H{9} -^ d6, (96) 



where 7 = 2(1 — H_). Because the right-hand side is bounded, this proves that the as- 
sumption 2(b) (in Theorem 1) is satisfied with L(h) = 1. 
Computing (Si o 8\ R)(s, t), we get 

(#o«y*R)( Sj t) 

A(6)(-\s - t - 2h\ 2H ^ + i\s - t - h\ 2H W - 6\s - t\ 2H ^ 

+ 4\s - t + h\ 2H ^ -\s-t+ 2h\ 2H W) dO. 



Thus 



{S'to5%R)(t,t) = 8 / k(9)\h\ 2H[e) {A-2 2H{e) )&e, 



^ °^ 2 R l {t,t) = 8 j\(e)\h\ 2 ^-^(4 - 2 2g W) c\9 (97) 



with 7 = 2 - 2H. 

Setting g (t) = (4 - 2 2 ^)J K , we have 



\h\ 



2-7 



■go® = 8 rAWl {ffW>m (4-2 2 ^))|/ 1 | 2 ^)-^d0. (98) 
Jo 



Therefore, the Lebesgue theorem and the fact that the right-hand side does not depend 
on t yield 

(#0 i) 



lim sup 

h->0+ h<t<l-h 



\h\ 



2-7 



9o{t) 



0. 



Hence the assumption 2(c) (in Theorem 1) is fulfilled. 

So if we apply Theorem 1 to Z, wc obtain (79). □ 

Proof of Proposition 11. We apply Theorem 2 to Z. We must show that the assump- 
tions 2 and 3 (in Theorem 2) are satisfied. For assumption 2 (in Theorem 2), we must 
compute the derivative gfr^r- As in the proof of Proposition 10, the Lebesgue theorem 
of differentiability under the symbol integral shows the existence of this derivative on 
]0, l] 2 \ {s = t} and yields the formula (96), 

= -4 / A(6)2H(6)(2H(9) - 1)(2H(0) - 2)(2H{6) - 3)|s - t\ 2( * H ^~^ dO, 
Jo 
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where 7 = 2(1 — If). Hence the assumption 2 (in Theorem 2) is satisfied with L(h) = 1 
and C(s, t) equal to the right-hand side of (96) (the continuity of Con{s = i] is obtained 
by applying Lebesgue theorem). 

Next we deal with the assumption 3 (in Theorem 2). We have proved (97): 



(ff 0<%i?)(t,t) 
1/^-7 

We set (7 = 0, uq = 1/2 and 



A(6)\h\ 2 ( H ^~^(4-2 2H W)d9. 



g (t) = 8(4 - 2 2 ^) / A(0)1 {HW =h } d0 = (4 - 2 2 ^) % . 
Jo 



We have 



- = 8 /" m\h\ 2(H(6) - m (4 - 2 2ff W)l 
l"l 7 Jo 



{_f/(e)> ; H+i/4} 



d0. 



Therefore, the Lebesgue theorem yields 



lim —= sup 

fi->0+ V/l h<t<l-h 



«o^i?)(M) 



N 2 ' 



= 



and assumption 3(d) (in Theorem 2) is fulfilled. 
Moreover, we have 



so 



$o5$R)(t + h,t)=4 ! / A(0)\h\ 2H W(4.2 2H ^ -3 2H(e *> -7)d9, 



{5>to8%R){t + h,t) 



\h\*-t 



- 4 [" A(6)\h\ 2 ^ H ^'^(4.2 2H ^ - 3 2ff W - 7)d0. 
Jo 



Consequently, the theorem of dominated convergence yields 

{6$o5%R){t + h,t) 



lim sup 

fi->0+ h<t<l-h 



\h\ 



2- 7 



0. 



where 



3(;) = i(2 2 ^+ 2 -3 2 ^-7)JH. 
Thus, the assumption 3 (in Theorem 2) is fulfilled. 



□ 



Proof of Corollary 12. The almost sure convergence is a straightforward consequence 
of (79). To get (82), we just apply the 8 method with the C 1 function: 



f(x,y) 



\og{y/x) 
21og2 ' 



□ 



Quadratic variations of Gaussian processes 



749 



Proof of Proposition 13. Wc apply Theorem 1 to Z to obtain a refinement of Propo- 
sition 10. We must show that Assumptions 2(b), 2(c) and 3 (in Theorem 1) arc satisfied. 
For assumption 2(b) (in Theorem 1), we use the same arguments as in Proposition 11 
and obtain the same function C(s,t) with 7 = 2 — 2H_ and L(h) = 1. 

Next we deal with assumption 2(c) (in Theorem 1). Assumption 2(c)(iii) (in Theorem 
1) is a straightforward consequence of the Lebesgue theorem. 

Moreover, we have proved (97): 

(<°g?)(M) = s£ A(9)\h\ 2 ^-^(A - 2 2H ^)d9. 

We set q= 1, ct\ = 1/2, V\ = \ and 

g (t) = 8(4 - 2 2 ^) r A(0)l mg)=m d6 = (4 - 2 2 ^), 
Jo 

9i(t) = 1, 

4>{h) =8 rA(0)(4-2 2 ^))|^| 2 W)-^l {ii < ffW <H + i/4 } d0. 
Jo 

We have 

v l^p.^ - .9oW - gi(t)4>(h) 

= 8 rA(0)(4-2 2 ^))|/ 1 | 2 (^)-^l m+1/4<ff(e)} d0 
./o 

,l =° + o(^), 

thanks to Lebesgue theorem. Because the right-hand side does not depend on t, assump- 
tion 1(c) (v) (in Theorem 1) is fulfilled. 
Moreover, we have 

= 8 rA(0)(4-2 2 ^))|/ l | 2 W )-Sl m<HW<ii+1/4} d0. 
Jo 

Therefore, for h enough small, 

</>(h) > Kh 2 ^ +1 ' A -^ = KVh, 

which yields that the assumption 2(c) (iii) (in Theorem 1) is satisfied too. The last in- 
equality yields that the assumption 3 (in Theorem 1) is fulfilled. 

Therefore, (83) is a consequence of Theorem 1 applied to Z. To prove (84), we ap- 
ply Theorem 2. We need to check that assumptions 3(d) and 3(e) (in Theorem 2) are 
fulfilled. Assumption 3(d) (in Theorem 2) is a straightforward consequence of previous 
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computations. For assumption 3(e) (in Theorem 2), we use the same arguments and the 
same function g(t) as in Proposition 11. □ 



Proof of Corollary 14. The almost sure convergence is a straightforward consequence 

n 2HK-i Vn{z) 

K(2n) 2HK - 1 V 2n (Z)) 



of (79). To prove (86), we apply the S method between the points L 2n ^ 2 HK-iy / Z y an< ^ 



{^tiy^k \) to the C 1 function: 



f{X ' y) ~ 21og2 ■ □ 



4.4. Proofs of the results in the non-1. a.s.s. case 



Proof of Proposition 15. Thanks to Proposition 9, we can compute the variogram of 
A^ f or all teR 2 : 

«(*) = § / C{l,H(e)) 2 \t lC os(d)+t 2 sm(9)\ 2HW d9. 



We use the polar coordinates and we parametrize tsR 2 \ {(0, 0)} by (p(t), a(t)): 

t\ = p(t) cosa(i), 
^2 = p(t) sina(i). 

We set v e (t) = e~ 2 —v(et),e > 0, and use the polar parametrization and the 7t-periodicity 
of the function H to obtain 

v e {t) = \ r +n/2 C(l,i7(0))VW 2ff(e) |cos(a(t)-0)| 2ff(e) e 2(ff W^d0. 

We assume that p(t) ^ and a(t) ^ 6* + Tt/2 (mod7t). The Laplace method (see 
Dieudonne [11], Theorem IV. 2. 5, page 125) yields 

vs(t) ^#^)^|cos(a(i) 6^, I W 



{-\oge)H»{6*y 
Therefore, 



lim J log I l)-L v (et)= C(1 ;~ )2 p(t) 2 H\cos(a(t) ~ 6*)^ ^ 



6 _o+ y °\e J e 2 " v ' 8 rw 1 v w " V H"(p.) 

With a refinement of the Laplace method, we can check that this is still true if p(t) = 
or a(t) = 0* +7t/2 (mod7t), which are the cases where the limit is vanishing. 
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Moreover, because A^ H > is a Gaussian processes with stationary increments and van- 
ishing limits at the origin, we have 

Vs, t G R d Cov(4 ff ) , A[ H) ) = 2(v(s) +v(t)- v(s -t)). 

We are able to conclude that the finite-dimensional laws of the process (87), with to = 
0, converge toward those of the process defined in (88). Because has stationary 

increments, we have the same result for all to G R d - □ 

Proof of Proposition 16. As in Proposition 10, we must check that the assumptions 
of Theorem 1 are satisfied. For assumption 2(b), in Proposition 10 we showed formula 
(96), 

= -4 / h(6)2H(0)(2H(6) - l)(2H(9) - 2)(2H(9) - 3)|s - t\ 2{H ^S- ] d6, 
Jo 



where 7 = 2(1 — H_). We set L(h) = J j^y L , where L is slowly varying, and write 

L(\ s -t\) ds 2 dt^ s,t > 

k(0)2H(0)(2H(6) - 1)(2H(0) - 2){2H{6) - 3) 



x y/\og(\s-t\)\s - t|aW»)-S dO. (99) 
If H ^ 5 , the Laplace method yields that the right-hand side goes to 



A(6*)2H(2H-l)(2H-2)(2H-3) (100) 



H"{0,) 



when \s — 1\ — > 0. A refinement of the Laplace method allows us to check that it is still 
true if H_ = h. This implies that the right-hand side of (99) is a continuous function of 
(s,t) on the set [0, l] 2 . Therefore, it is bounded and the assumption 2(b) (in Theorem 1) 
is satisfied. 

For assumption 2(c) (in Theorem 1), we have seen in (97) that 

vjhjMM =B> j\{6)\h\^ H ^S-\A - 2 2H ^)de 

with 7 = 2(1 — IP). We want to study the asymptotic behavior of the preceding integral, 
so we denote by 1(h) the right-hand side of (97). The Laplace method yields 



h^0+ „ k rn w„ ~ 2 H\ / 27t 



1(h) ~ 8A(0»)(4-2^). 



■\ogh)2H"(G*)' 
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Therefore, we obtain 

h Hm y/(-l gh)I(h) = 8A(5»)(4 - * K )sJh^- 

As in the Proof of Proposition 10, because y/ (— log h) 1(h) does not depend on the variable 
t, we have 



lim sup 

h^0+ h<t<l-h 



(6*o6*R)(t,t) 

flow 



L(h)\h 



2-7 



:(). 



where go(t) = (4— 2 2 —)Gg r . Therefore, the assumption 2(c) (in Theorem 1) is fulfilled 
with g = 0. □ 

Proof of Proposition 17. We apply Theorem 1. The difference from the proof of Propo- 
sition 16 comes from assumption 2(c) (in Theorem 1). We must compute an asymptotic 
expansion of the expression, which is a consequence of (97), 



(5<to5 h 2 R)(t,t) 



1-271 

■■4y/=to£h / A(0)|ft| a ™- m >(4-2 2jr W)d0 
Jo 



m\h\*-i 

with 7 = 2(1 — H). The 7t-periodicity of the functions H and A yields 

(ff o ffy t) = 8 yn^ f 9 ' +n A(9)(4 - 22*W) e a™-«Oiogh d 0. 

L(/l)|/l| 2 "T 7 e-> 

Here we cut this integral into two parts (we integrate on [0*,6* + n/2] and on [6* + 
n/2, 8* + n]) and use Theorem II. 1.1 from Wong [17] on the extended Laplace method. 
We obtain 

(6*o6*R)(t,t) /^p+ ^ r(^ + l)/2)a, / 1 \ 
L(h)\h\ 2 ~y ^ (-logh)*/a + {(-logh)"/ 2 ]' [ 1 

where the coefficients <7j can be expressed in terms of Hi and Aj. 

Because these quantities do not depend on the variable t, the assumption 2(c)(iii) of 
Theorem 1 is fulfilled with 

g l (t) = 16r( l -^-)a l , £, = 1/2, <f>(h) 



Proof of Proposition 18. We apply Theorem 2 to Z. We must show that assumptions 
2 and 3 (in Theorem 2) are satisfied. As in Proposition 16, we can check that assumption 2 
(in Theorem 2) is fulfilled with C(s,t) equal to the right-hand side of (99). 
For assumption 3 (in Theorem 2), we have proved (97): 

= 8 J* m \ h \2(H(8)-H) {4 _ 2 2H ^)d9. 
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We set q = 0, oi\ = 1/2, V\ = 1 and 

g (t) = (4-2 2 ^)Gg„ 
Si(*) = l, 

0(h) = Sv^loiTI f A(0)(4 - 2 2ff ( e ))|^| 2 ( ff ( e )-^l {ff(e)<ii+1/4} dfl 
Jo 

-(4-2 2 ^)G e ,. 

We have 

= 8V r I^rAW(4-2 2 ^))|/ l | 2 ^)-^l {£+1/4<H(e)} d^ 



o 





thanks to the theorem of dominated convergence. Because the right-hand side does not 
depend on t, assumption 3(d) (in Theorem 2) is fulfilled. 
Moreover, we have 

(<ff o 5%R)(t + h,t)=4 f k{9) \h\ 2H ^ (A.2 2H ^ - 3 2ff W - 7) d6 
Jo 

and so 

^VSa- M) = 4 v /T ^ f A(fl)|h| J ( fl W-S(42 aif W - 3 2H ^ - 1)60. 
L(h)\h\ 7 Jo 

Consequently, the Laplace method and the theorem of dominated convergence yield 
(5>?o6$R)(t + h,t) Ak(a . tAtt2H o2H „ r^T 



lim sup 

/i->0+ h<t<l-h 



4A(6>*)(4.2^--3 2 ^-7) 



= 0. 



L(/i)|/i| 2 " 7 v ' 'V 

Thus, the assumption 3(c) (in Theorem 2) is fulfilled with 

^2^-3^-7 
ffW = g e * ■ D 
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